Effect of quark masses on the QCD presssure in a strong magnetic background 
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We compute the two- loop contribution to the QCD pressure in a strong magnetic background, 
for arbitrary quark masses. We show that, for very large fields, the chiral limit is trivial. 



Large magnetic fields can be created not only in the core of magnetars [l[ but also in current experiments at 
BNL/RHIC and CERN/LHC involving non-central heavy ion collisions. The fields created in these collisions are pos- 
sibly the largest magnetic fields produced since theprimordial electroweak transition, reaching values B ~ 10 19 Gauss 
{eB ~ 6 m 2 ) for peripheral collisions at RHIC 2], and even much higher at the LHC thanks to the fluctuations 
in the distribution of protons inside the nuclei [3|. Such intense magnetic fields may dramatically affect the phases 
of strongly interacting matter, as is the case in more ordinary circumstances [4(. The mapping of the QCD phase 
diagram in the T — eB plane is still in its infancy (see e.g. [5| and references therein). There are clear indications 
that sufficiently large magnetic fields do modify the nature and behavior of the chiral and the deconfinement phase 
transitions [tML9j. New phases are also predicted [20l-[23|, and it has even been suggested that the vacuum may turn 
into a superconducting medium via p- meson condensation [24j . 

While most of the analyses so far have relied on effective models, or calculations in large N c limit of QCD [25], the 
first results from lattice QCD have been obtained recently [26l-[29j. This opens a new channel for comparison between 
analytical or semi-analytical techniques and numerical non-perturbative approaches. In this perspective, we note that 
the recently established discrepancy between different lattice QCD results (for large [l?} and physical values [H, l29j 
of quark masses) is most likely related to quark mass effects. It is the purpose of this paper to analyze the possible 
competition between mass and magnetic-field corrections to the QCD pressure. More specifically, we compute the 
two- loop correction to the QCD pressure in a magnetic background field and for arbitrary quark masses. We indeed 
find a significant competition between the effects of quark masses and those of the magnetic background. In particular, 
for extremely intense magnetic fields, we show that the two- loop contribution to the pressure is trivial in the chiral 
limit. 

We shall assume in our calculation a constant and uniform Abelian magnetic background, whose strength is large 
enough to produce interesting effects, i.e. eB > m 2 . We also consider the temperature to be large enough that 
perturbation theory can be applied to the calculation of the pressure. Nevertheless, the effect of the strong background 
magnetic field must be treated non-perturbatively: this is achieved by using the propagator that was obtained long ago 
by Schwinger (30[, and that can be cast in a convenient form using Landau levels, as shown in Ref. (31] (see also Refs. 
[3214341]). Since we restrict our analysis to the case of very intense magnetic fields, the summation over the Landau 
levels is rapidly convergent, and the leading correction to the pressure is obtained from the Lowest Landau Level (we 
shall refer to such calculations as the lowest Landau level (LLL) approximation). The corresponding propagator for 
a fermion of a given flavor / and (absolute) electric charge g/, in the presence of the classical field A c \ = (0, A) (with 
V x A = B = Bz) reads: 
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where we have used a compact notation for the transverse (pt = (^1,^2), 7t = (7\7 2 )) an d longitudinal (pl = 
(Po 5 P3) 5 7l = (7° 7 7 3 )) quantities. This is equivalent to the result used in Ref. [33|, obtained by constructing the 
projectors on the different Landau levels from the exact solution of the Dirac equation [35|. A peculiarity of the LLL 
approximation should be noted: the propagator is a 4 x 4 matrix, but it describes only two physical propagating 
modes. Two eigenvalues of Sq LL indeed go to zero in the vicinity of the lowest Landau level pole (pfj = m 2 f 
This complicates in particular the computation of the free pressure in the LLL approximation [l8|, [34| . 

The thermodynamic potential of QCD, up to two loops, is obtained form the standard diagrammatic expansion: 



2 





[diagrams with counterterms] + 0(3 loops), 



(2) 



where full lines are fermions, dressed by the magnetic field, curly lines are gluons and dashed lines represent ghosts 
(whose role is essentially to cancel the contribution of spurious degrees of freedom in the gluonic pressure). The 
calculation is carried out in Feynman gauge. 

The gluonic part is equivalent to the usual hot perturbative QCD result and is, therefore, well-known (36^ : 
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The one-loop contribution to the fermionic pressure has been considered in different contexts (usually, in effective 
field theories 0-EI Qj| Gil HtJ) an d computed from the direct knowledge of the Landau levels E 2 (n,ps) = p\ + mj + 



2qfBn and their degeneracies qfB/(2ir) for n = and qfB/ir for n = 1,2, • 



The final exact result reads (see Ref. 



for discussions on the subtraction procedure) 
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where fif is the quark chemical potential (associated to baryon number conservation). In the limit of large magnetic 
field (i.e. Xf = m^/(2qfB) — )> 0), reduces to the LLL contribution 
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The exchange diagram corresponds to the first nontrivial contribution. In terms of the propagators in coordinate 
space, this diagram is given by 
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where A a are Gell-Mann matrices, with A a A a = — l)/2, the trace Tr acts over Dirac indices and the 4-momentum 
is given in terms of the Matsubara frequencies (ouf = 2lnT) and of the 3-momentum k as: K = (k° = iujf , k) . 

Notice that if we assume translational invariance and the free Dirac propagator for the fermions (with P = 
(p° = iuo^ + •> p) an d u n = (2n + l)7rT), this expression reduces to the usual one [36|. In the presence of a 
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uniform and constant magnetic background (B = Bz), however, the fermion propagator becomes dependent on x and 
y in a nontrivial way due to the Schwinger phase, as discussed previously. 

However, a detailed analysis of this diagram shows that it can be cast in the following neat form [34 
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where (y c = N C (N 2 C - l)/2) 

I , TL 1 , 7T<2 

and \zl = (ic^P, fc^), Pl = (iw^iPz), Pl — {iuj^q z ). This expression realizes concretely the intuitive expectation that 
the nontrivial dynamics in a strong magnetic field is effectively (1 + 1) dimensional 1 . Since the gluons do not couple 
directly to the magnetic field, their dispersion relation maintains its (3 + l)-dimensional character (uj 2 = k 2 + k 2 + &f), 
which effectively results in a "massive" gluon (m 2 = k 2 -\- k 2 ) in the reduced (1 + l)-dimensional diagram. In 
the end the exchange contribution to the QCD pressure in the lowest-Landau level approximation for the fermion 
propagation is essentially an average over the effective gluon transverse mass m 2 = k 2 + k 2 of the exchange diagram 
in (1 + l)-dimensions with the Gaussian weight (qfB/27r)ex.p[—m 2 z /2qfB]. 

Apart from the tensorial structure, this diagram corresponds to the two-dimensional version of the exchange diagram 
of a Yukawa theory with both massive fermions and bosons, which was computed originally in Ref. [38j (cf. also 
[39|| ) . In the present (1 + l)-dimensional context, however, one expects renormalization to be trivial. It can be shown 
indeed that the usual quark self-energy counterterms vanish in dimensional regularization, while the high momenta 
in the gluon lines are fully tamed by the Gaussian weight (qfB/2n)exip[—rri 2 ; /2qfB]. Our concern here will be thus 
directed towards the IR domain, which proves to be subtle as discussed for instance in Ref. [331 ]. 

Since the result in Eq. ([7]) is (superficially) proportional to m 2 (see Eq. [8j), it vanishes in the chiral limit rrif — >• 0. 
Simple power counting suggests that this result is not affected by the remaining integrations over the longitudinal 
momenta, which are are only logarithmically divergent when rrif — >• 0. This is confirmed by a numerical analysis of 
the complete integral. This is calculated by performing first the sums over the Matsubara frequencies, via standard 
contour integration. We obtain then, for one massive flavor [34| 
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where n B and rip are the Bose-Einstein and Fermi-Dirac distributions, respectively, uj = ^k 2 + , and 
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The transverse motion is not affected by the temperature, and is dictated by the lowest Landau level. Increasing the magnetic field 
shrinks the individual Landau orbits as r ~ 1/ ^qfB. So, the orbital motion of quarks in the plane transverse to the magnetic field 
becomes more and more constrained as B grows. In the limit of very large fields, the original helicoidal (tubular-like) paths become 
essentially straight lines parallel to the field direction. Of course, the longitudinal motion is affected by the heat bath, so that the 
tubular structures become "blurred" (noisy). However, if the magnetic field is larger than the temperature, this effect will be minor and 
the quark motion will be essentially one dimensional. 



The integrals over can be written in radial coordinates and put in a more convenient form by the change of 
variables x = |kj_|/ \2qfB). The Gaussian weight becomes then a representation of the Dirac delta function for very 

C V a/ 71 "/^ ~^ d(x)), so that the UV limit is completely under control. By 



large magnetic fields (for b = 2qfB, e 

performing the integrals numerically, one finds that the overall mj factor controls indeed the IR sector: without this 
factor, the integrals would diverge in what seems to be a logarithmic fashion (as can be also seen from semi- analytic 
stronger approximations in the limit of large fields). This confirms the result anticipated in Eq. (j8j): the exchange 
contribution to the QCD pressure for very high magnetic fields vanishes in the chiral limit. 

It is clear, then, that the quark masses play an important role in the perturbative calculation of thermal QCD 
under an external magnetic background, competing with the field. In fact, although the effects from quark masses 
on the thermodynamics of QCD have been greatly overlooked for many years, it has been shown that they can 
bring significant corrections to the perturbative pressure at finite density [40j |. with co nseq uences to the structure of 
compact stars, and at finite temperature [4l|, affecting quark mass thresholds (see also |42| ) . In order to make direct 
comparisons to our results at large magnetic fields, we also present here an explicit analysis of the mass-dependence of 
the 0(g 2 ) thermal QCD pressure. The derivation follows again the same steps as the analogous computation for the 
massive Yukawa theory [38] in dimension 4, with the simplification that gluons are massless. For one-flavor QCD, after 
renormalization (in the MS scheme), one obtains for the exchange contribution to the thermodynamical potential: 
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FIG. 1: Exchange contribution to the pressure, p^ xch ; as a function of the quark mass for qf — 2/3, T = 100 MeV and 
eB/m% = (lower), 100, 200 and 300 (upper). The strong coupling is fixed at a s = g 2 /47T = 0.3 and the renormalization scale 
for the zero magnetic field curve is set to A = 800 MeV. 

In Figure [1] we display the behavior of the contribution from the exchange diagram to the pressure as a function of 
the quark mass for qf = 2/3, T = 100 MeV and a few values of eB in units of the pion mass, as well as for vanishing 
eB. The picture clearly shows that the exchange contribution vanishes smoothly in the chiral limit. As the quark mass 
is increased, the exchange pressure becomes nonzero but remains always between one and two orders of magnitude 
below the leading (free) contribution 2 . This fact might indicate an improved convergence of the perturbative series 
in the quark sector 3 in this limit of high temperatures and extremely large magnetic fields. This is further supported 



2 A thorough comparison between the free pressure and the 0(a s ) contribution as well as thermal effects will be presented in a longer 
publication. 

3 Notice that the purely gluonic sector remains unaffected, so that it is obviously still plagued by infrared divergences as its B = 
analog. Nevertheless, for extremely large fields, the contribution of pure-gluon diagrams to the pressure (which are ^-independent) can 
be neglected as a subleading correction and the QCD pressure will be dominated by the quark sector. 
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by the absence of any explicit dependence of the pressure on the renormalization scale. Higher-loop computations or 
a direct comparison with lattice QCD data for the pressure could confirm this in the near future. 

The comparison with the behavior of the pressure at zero magnetic fields shows that the exchange contribution to the 
pressure changes its sign when the QCD medium is exposed to a very intense magnetic background, \fe~B ^> T. This 
qualitative change is directly related to the effective dimensional reduction at large magnetic field and the consequent 
modification of the Dirac traces (the vanishing of the exchange pressure in the chiral limit can be understood as 
resulting from a conflict arising, in one dimension, from the fact that the exchange diagram couples right and left 
movers, while the vertices conserve helicity). The nonzero effective mass acquired by the gluon in the dimensionally- 
reduced diagram also contributes to prevent cancellations that usually happen at zero magnetic fields. 

We would like to conclude with an intriguing observation. It is well-known that the gyromagnetic factor of Dirac 
fermions, g m , is not exactly 2, but that it receives radiative corrections from QED (also from QCD, but these are 
subleading). The deviation, although small (g — 2 ~ 10 -3 ) may produce sizable corrections to the pressure. Indeed, 
it affects the energy of the lowest Landau level, effectively turning the mass m into m^ff = m 2 + (g — 2)eB. For 
eB ~ m 2 , the correction is in the MeV range; it may thus compete with m and cannot be ignored in a quantitative 
study. 4 
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